June 1996 
DESY 96-100 
hep-ph/9606242 

The Coefficient 773 of the |AS| = 2-Hamiltonian in the Next-To-Leading 
Order 

Stefan Herrlich a * 

a DESY-IfH, Platanenallee 6, D-15738 Zeuthen, Germany 

I present the calculation of the QCD short distance coefficient 773 of the | AS| = 2-hamiltonian in the next-to- 
leading order (NLO) of renormalization group improved perturbation theory. It involves the two-loop mixing of 
bilocal structures composed of two | AS| = 1 operators into | AS| = 2 operators. The next-to- leading order corrections 
enhance 773 by 27% to 



thereby affecting the phenomenology of the CP-parameter 6k sizeably. 773 depends on the physical input parameters 
m t , m c and A^g- only weakly. The quoted error stems from factorization scale dependences, which have reduced 
compared to the old leading log result. We further discuss some field theoretical aspects of the calculation such 
as the renormalization group equation for Green's functions with two operator insertions and the renormalization 
scheme dependence caused by the presence of evanescent operators. This article is based on work done in 
collaboration with U. Nierste. 



1. Introduction 

The effective low-energy hamiltonian inducing 
the I AS I = 2-transition reads: 



H 



AS 1=2 



r 2 

+2X c X tV * 3 S(x^x* t ) 
+h.c. 



2X c XtH ci 
X 2 cV lS(x* c ) + X 2 t ^S(xt) 



(1) 



Here Gf denotes Fermi's constant, Mw is the W 
boson mass, Xj = VjdV* s ,j = c,t comprises the 

CKM-factors, and Qs2 is the local dimension-six 
|AS| = 2 four-quark operator 



/S2 



= [s7m (1 - 75) d] ■ [sjf, (1 - 75) d) 



(2) 



m* /Mjy, q — c,t encode the running 
m q {m q ) in the MS scheme. In 
writing (0) the GIM mechanism A u + X c + Xt = 
has been used to eliminate A u . Further we have set 



The x* 
quark masses m* 
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m u = 0. The Inami-Lim functions S(x), S(x,y) 
describe the | AS | = 2-transition amplitude in the 
absence of QCD. They read: 

"19 1 3 1 
4 



S(x t ) 



S(x c ) 
S(x c ,x t ) 

with 

F(x t ) = 



1 - x t 
,2 



41 -x t 

3 

lnx t , 



2(1- x t ) 



x c + 0{xi), 
-x c \nx c + x c F(x t ) 
+0{x 2 c ]nx c ), 



x 2 t 



8x t +4 



In a; t 



XL 



(3a) 
(3b) 

(3c) 
(4) 



4(1 -x t ) 2 ' Axt-1 

In © and © we have only kept terms which are 
larger than those of order (m s m c ) /M$y neglected 
by setting the external momenta to zero. 

In (Q) the short- distance QCD corrections are 
comprised in the coefficients and 773 with 

their explicit dependence on the renormalization 
scale /i factored out in the function b(fi). In ab- 
sence of QCD corrections r)*b(fi) = 1. 
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Table |l| summarizes the logarithms summed by 
the forthcoming renormalization group (RG) evo- 
lution from My/ down to m c in the different or- 
ders. 

Table 1 

Logarithms summed by the RG evolution from 
Mw down to m c for the three terms in (|IJ), n = 
0, 1, 2, . . . The last line shows the order in which 
the dependence on m t enters. 



Order 


H c 


H l 


H ct 


LO 


(a s \nx c ) n 


(a s \nx c ) n 


(a s lnx c )"lna; c 


NLO 


a s (a s \nx c ) n 


a s (a s \nx c ) n 


(a s lna; c ) n 


mt 


none 


in LO 


in NLO 



The first complete determination of the coeffi- 
cients i = 1,2,3 in the leading order (LO) is 
due to Gilman and Wise [Q. However, the LO 
expressions have several conceptual drawbacks: 

i) The fundamental QCD scale parameter 
Aj^g is not well-defined in the LO. 

ii) The quark mass dependence of the 7/i's is 
not correctly reproduced by the LO expres- 
sions. Especially the m^-dependent terms 
in 773 • S(x*, xX) already belong to the NLO. 

iii) Similarly the question of the definition of 
the quark masses (i.e. the renormalization 
scheme and scale) to be used in (]l|) is a next- 
to-leading order issue. 

iv) The LO results for 7/1 and 773 show a large 
dependence on the factorization scales, at 
which one integrates out heavy particles. 
In the NLO these uncertainties are reduced 
considerably. 

v) One must go to the NLO to judge whether 
perturbation theory works, i.e. whether the 
radiative corrections are small. After all the 
corrections can be sizeable. 

To overcome the limitations listed above one has 
to go to the next-to- leading order (NLO). This 
programme has been started with the calculation 
of r\\ by Buras, Jamin and Weisz Q . Then Nierste 
and I have derived the NLO expressions for r\\ [|| 



and 773 Q. This article will essentially deal with 
the term H ct of (|l|) and is based on Q . 

2. The NLO calculation of 7/3 above the 
charm threshold 

Here we will shortly describe how the large loga- 
rithm In x c present in ( |3c| ) is summed to all orders 
in perturbation theory. This is done in two steps: 
First one sets up an effective lagrangian C cS in 
which the W boson and the top quark are removed 
as dynamic degrees of freedom. In /Cj,^ 8 ' -2 the 
I AS I = 1 and | AS| = 2 transitions are described by 
local four-quark operators, which are multiplied 

by Wilson coefficients. The general structure of 
£|AS|-2 rea( j s: 

£^ S ' 2 = 7=f VcKM ^ CkQk 

V 2 k 

-^fVcKM^QQi. (5) 

Here the Vckm denote products of CKM elements. 
The Q k , Qi represent local |AS| = 1 and |AS| = 2 
operators and the C'k, Ci are the corresponding 
Wilson coefficient functions. The |AS| = 1 part 
of (||) contributes to |AS| = 2 transitions via dia- 
grams of the type displayed in Fig. [I]. The Wil- 




Figure 1. The diagrams Do and Po in the effec- 
tive five- and four-quark theory. The light crosses 
denote insertions of local |AS| = 1 current-current 
operators, the grey ones insertion of local | AS | = 1 
four-quark penguin operators. 

son coefficients are fixed at a factorization scale 
Htw = 0(Mw,Tnt) by the requirement that the 
I AS I = 1 and | AS | = 2 Green's functions derived 
from (|^) are equal to the same quantities cal- 
culated with the full SM lagrangian. The ef- 
fect of this procedure is that the lnx c present 
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in (pq ) is split as \a.{jjL t w I ^Iw) + hi(m c //itiv). 
Here the former term resides in the Wilson co- 
efficients while the latter is contained in the ma- 
trix elements of the local operators. The choice 
fj-tw = 0(Mw, m t ) ensures that the Wilson coef- 
ficients do not contain large logarithms and there- 
fore can be reliably calculated in ordinary pertur- 
bation theory. 

The second step is the RG evolution of the 
Wilson coefficients from the scale faw down to 
fj, c = 0(m c ) which sums ln(fj, c / fx t w) -0 

2.1. The Operator Basis 

Let us now construct C^ S ^ in (||) as far it is 
needed for the calculation of 773. Since the pres- 
ence of |AS| = 2 terms in (||) does not affect the 
|AS| = 1 part of Z^ 8 ' -2 , we can simply take the 
latter from They consist of the following 

set of operators: 



reads 



5 2 M 2e 



QS2 = 



S 2 M 2e 



s^Ld-s^Ld, (7) 



which follows from power counting and the ab- 
sence of any non-zero mass parameter apart from 
m c . The inverse powers of g are introduced for 
later convenience as in |Q. One may arbitrarily 
shift such factors from the Wilson coefficient into 
the definition of the operator. The factor yT 2e 
stems from ^barc = Z g gfi e and the fact that 



Q1 1 


= (sJn Lk ) 


■ (l^Ld) 


1, 




(6a) 


Q2 1 


= (sj^Lk) 


■ (l^Ld) 


1, 




(6b) 


Q 3 


= (sj^Ld) ■ 


E 

q—d,u,s : ... 


frfLq) ■ 


1, 


(6c) 


Qi 


= (sj^Ld) ■ 


E 

q—d,u,s,... 


frfLq) ■ 


1, 


(6d) 


Qs 


= (sj^Ld) ■ 


E 

q—d,u,s,... 


(jtfRq) 


• i, 


(6e) 


Qe 


= (sj^Ld) ■ 


E 




• i. 


(6f) 



q— d,u,s,... 

Here the Q\ l , i = 1,2, k,l = u,c represent the 
|AS| = 1 current-current operators, the Qi, i = 
3, . . . , 6 the QCD-penguin operators. The sum 
in (|6c|)-(|6|) runs over all active flavours. Further 
L,R = (1 =F 75) and 1, 1 denote color singlet and 
anti-singlet, i.e. Qf = (sij^Lkj) ■ (Jjj^Ld^j with 
i,j being color indices. 

Further we need the |AS| =2 operators present 
in (||). Consider first the diagram Do of Fig. [j] 
with two internal charm quarks and zero exter- 
nal momenta. The only physical dimension-eight 
I AS j = 2 operator required to absorb its divergence 



^For simplicity we ignore the intermediate scale fib 
0(m,i,) at which the bottom quark gets integrated out. 



-ibarc 

h 



^barc 



[sj^Ld ■ sj^Ld] 



(8) 



must be independent of /i^] Any other dimension- 
eight I AS I = 2 operator contains one or two powers 
of m c less than Qj and derivatives and/or gluon 
fields instead. Their on-shell matrix elements are 
suppressed by powers of m s /m c with respect to 
those of Q-j, so that they do not contribute to the 
coefficient of the leading dimension-six operator 
below the charm threshold (cf. ([j])). Likewise they 
cannot mix with Q7 under renormalization. 
Therefore our operator basis consists of Qi l 2 , 



2 ) Here and in the following the superscript "bare" denotes 
unrenormalizcd operators, while renormalized ones do not 
carry an additional superscript. 
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Qz 6 and Q7^]and C cS S is found as 



H,bare 



J=l k,l—u,c 



bare 



i=3 



G7 



2 6 



Q^ arc 



-counterterms proportional to (12) 
unphysical operators. 



3 'ln addition to these physical operators we have to take 
into account several evanescent operators. This class of 
operators appears quite naturally when one has to deal 
with the renormalization of operators containing rnore 
than one fermion line in dimensional regularization [H. 8|f| . 
To illustrate some of the findings of |J we have defined the 
evanescent operators with some arbitrary coefficients ai, 
a 2 , 01, bi: 



ExlQj] 
Bi[Q 7 ] 

E 2 [Q 7 ] 



- (4 + die) 7m L ® 7 m l] K xj , j = 1, . . .4, (9a) 

- (16 + a 2 e) 7 M i? ® Kij, j = 5, 6, (9b) 



— (4 + aie)7 M L(g)7 fl Ll #12, 



[7n7^7i77CT7ri ® 7 T 7 CT 7''7"7 M L 
- [(4 + aic) 2 + 6ie] 7 M L <S> -y^L] K 22 . 
vitli the color factors 



K 12 = K13 = K 15 
Ku = K XA = K 16 



|1 



W 1 ^ 



(9c) 



(9d) 



(10) 



4Af 



Apart from places where it is indicated we will always state 
the results corresponding to 



-16, 



(11) 



in order to comply with the standard choice used in Q 
|^,^| [t|] . Since NLO anomalous dimensions and matching 
corrections of physical operators do not depend on bi , we 
do not give a numerical value. Likewise we do not need 
the value of the colour factor K 22 . 



2.2. The Effective Lagrangian C^ S]r ^ at the 
Scale [itw 

Besides the operators we need to know their 
corresponding Wilson coefficient functions at the 
initial scale \itw ■ The ones for the |AS| = 1 case, 
Ci, i = 1, . . . , 6 can be taken from |3|-fj]]. It is 
important to note that only C2 starts at 0(a°), 
the others at 0(a\). This means that the NLO 
matching can be done solely with the diagram Do 
of Fig. Q with two insertions of Qi. One easily 
finds 





4tt ^ n M w 
(x t (fitw)) + 2) 



in LO 



(13) 



NLO 



where F(xt) is the top dependent part of S{x c , Xt) 
defined in (Ji|). The factor a s originates from the 
special definition of Q7 in (Q) . Note how the large 
logarithm In x c in ([k]) is split between the Wilson 
coefficient C7 and the matrix element. The NLO 
result in (JlJ) is specific to the NDR scheme with 

© 4) 



I AS I— 2 

2.3. Evolving down C cS from [i t w to [i c 
Next we have to evolve down the Wilson coef- 
ficients present in £j,g S ' -2 ( fl2|) from ^tw to \i c . 
For the | AS| = 1 functions Ci, i = 1, . . . , 6 this is 
achieved by standard methods ||. To calculate 
the running of the |AS| —2 coefficient C7 we need 
to derive and solve the corresponding RG equa- 
tion. From Ms^^iff 3 ' -2 = in ( [l2|) one finds 

M^-^7 Cu)=CV (Ai)777 

2 6 

+EE ft M ft, w^'j (1 4 ) 

k=l k' = l 



with the anomalous dimension tensor^ 
-in 



«s~(0) 
7fe«,7= T-7fc„,7 



4 ) The Wilson coefficient Cj depends on ai 



■In 



-8 In 



M w 



+ 4F(x t (mw)) - (6 + ai) 



5 ' "Anomalous dimension tensor" is clearly a misnomer 
and only used to distinguish -fkn,7 from ordinary anoma- 
lous dimension (square) matrices. 
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^ ^ hkk'$nn' + Skk'lnn'} 



Z k'n',7 Z ^ 



k>=X n' = l 



z 



77- 



(15) 



It is possible to solve the inhomogeneous equa- 
tion ([l4]) directly but this turns out to be inconve- 
nient for practical purposes. Instead we may com- 
bine the evolution equations of the |AS| = 1 and 
I AS I = 2 Wilson coefficients into a single matrix 
equation. This is made possible because the GIM 
mechanism ensures that at least one of the two 
I AS I = 1 operator insertions in diagrams like the 
ones displayed in Fig. [I] is of the current-current 
type. The current-current part of the |AS| = 1 
mixing matrix, i.e. the entries related to Q±, Q2, 
can be diagonalized exactly using the basis 



}f = \ {Q k 2 l ± Qf ) . 



(16) 



Then (|l4j) together with the RG equation of the 
|AS| = 1 Wilson coefficients splits into two inde- 
pendent inhomogeneous RG equations 

V^-Cf (/x) = ynC? (^)+7 ±fc ,7C ± (fi) C k (ji) .(17) 

Here the decomposition of CV^tiy) hito C^(ntw) 
is completely arbitrary provided one satisfies 

C 7 {Htw) = C+ (n tw ) + Of {ntw) ■ (18) 
This decomposition is then automatically pre- 
served at any renormalization scale. 

Each of the two equations in ( |l7j ) may be writ- 
ten as a 7x7 matrix equation, which may be 
solved by standard methods. We can even do 
better and collapse the two resulting 7x7 matrix 
equations into one 8x8 matrix equation: 
d 



a/j, 
with 



D 



(19) 



We now need to know the elements of the 
anomalous dimension tensor 7±i,7> i — 1, . . . , 6. 
They are obtained from the renormalization con- 
stants using the definition in (|l^) and expanding 
the quantities in there in powers of a s and l/e. 
Here it is important to include the finite renor- 
malization terms needed for the correct treatment 
of the evanescent operators (^|). To calculate the 
LO term t£ 7 one needs to know the l/e parts 
of the one-loop diagrams displayed in Fig. |l|, the 
NLO part 7^7 requires the evaluation of a set of 



two-loop graphs. We find: 
/ -16 \ 



~(0) 

r+7 = 



\ 



'+,7 



-32 
-16 
32 

16 / 

/ -212 \ 

-28 
456 



~(o) 
T-7 



V 



1064 
3 

832 
3 



V 



8 \ 

16 


-16 

/ 

276 \ 

-92 

520 
-216 

1288 
3 

J 



(21a) 



• (21b) 



As usual the NLO anomalous dimension tensor 
depends on the renormalization scheme. The re- 



sult (21b) corresponds to the NDR scheme with 



the definition of the evanescent operators corre- 
sponding to (|ll|).[] 

3. The NLO calculation of 7/3 below the 
charm threshold 



7 = 



~T 
7+,7 

~T 

\ ?-,7 





777 - 7+ 




(20a) 



T±,7 : 



777 - 7- / 
(7±i,7i 7±2,7, 7±3>7, 7^7)^5,7, 7±ei7) (20b) 
/ (7( M ) \ 



D{jjl) = 



C+(m)/C+(/x) 

y c 7 -(m)/c_m J 



(20c) 



°)ln our two-loop calculation we have kept 01,02,11 and 
61 in (few arbitrary yielding 



~+.7 



188 74 _ 
"I T a l 

100 M ai 



3 3 

1816 
3 

680 



HI 



3 

1576 
3 

1664 
3 



80 „ 
— ai 

-fai 

28 „ 



aA 



7"2 



130 - 
3 

ff ai 
f#ai 



3 "' 1 
- — ai 

28 - 
" 3 



ai/ 



(22a) 
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AS 1—2 

3.1. The Effective Lagrangian C s at the 
Scale /i c = 0(m c ) 
After integrating out the charm quark all de- 
pendence on m c belongs to the Wilson coefficients. 
This implies that the term involving Q-j in ( |l2| ) 
has to disappear from the effective lagrangian, be- 
cause Qr contains m c in its definition ([?]). Fur- 
ther the |AS| = 1 operators are neglected in the 
new effective lagrangian, because the matrix ele- 
ments of double insertions of these operators are 
at most proportional to m 2 rather than m\. We 
have already neglected such terms in all preceding 
steps. 

Therefore the new effective lagrangian to de- 
scribe the physics below /i c reads: 



C 



|AS|=2 
off 



G 2 F r 



16tt 2 



2 Cs2 (m) + >H C si (M) 



+x c x t df (m) 



Z S2 (M) 1 



Sbarc 
'S2 ■ 



(23) 



This lagrangian already resembles — _ffl As l =2 in- 
troduced in ([!]) . For the matching we have to set 
the Green's function derived from (|l^) and the 
one derived from ( p3| ) equal at the scale /i = /i c . 

Let us start with the matching of CV in the LO: 
since the definition of Qi (0) contains the factor a s 
with respect to Qs2 we develop an explicit inverse 
power of a s for the Wilson coefficient: 



C*S2 ^ (Mc) — 



ml (/i c ) 4tt 



-C 7 0*c) in LO{24) 



2 a s (fi c ) 

Diagrams containing double insertions are of order 
a® and therefore start contributing to Cg^ in the 
NLO. The NLO version of can be written 



Cs2 ' (Mc) 



TO 2 (/i c ) 



S2 

1 4tt 



2 a s (Mc 



;C 7 (/i c ) 



*=+>- 



X] f «.S2 (/Uc) C-i (Mc) Cj (Mc) 

i=i 



(25) 



124 . 22 

~r+T ( 

-12 +6ai 

1496 
3 1 3 1 

-120 +8ai 



-oi 



+#ai 



A 



-4ai 



-4ai 



1160 
3 

-128 



16 „ 

-o-ai- 



-4ai 



"3 a2 + T ai 
-4a 2 -4ai/ 



(22b) 



The coefficients ri^si (a*c) in (l^) ar e given by 
the finite parts of the diagrams in Fig. [l| We find: 



Th^-l 



, m r (u r ) 



I In 



Mc 

m c {Hc) 

fJ-c 



nj for j =1,2, 
7ij forj =3,<26) 



nj for j =5, 6, 



where the r,-, 's denote the colour factors 



l±N 



T±l 

T+2 



T±3 = T± 5 



r +4 = 7+6 = 1, 
T-2 = r_ 4 = T_ 6 = 0. 



(27) 



Note that r^ga f° r J = 1,2 depends on the defini- 
tion of the evanescent operator £i [Q7] . As usual 
( |26| ) only holds in the NDR scheme. 

3.2. Evolving £^ S_2 below /i c 

The RG running in the effective three quark 
theory is particularly simple because there is only 
one operator left: <5s2- The running of the corre- 
sponding Wilson coefficient function Cgj can be 
taken from 

r/3 can then be determined from the identifica- 
tion 



/. I AS |=2 _ _ rr|AS|=2 
^cff — n 



(28) 



with C l £f^ from @ and ffl As l= 2 from (@). The 
result of this can be found in pL section 5] . 



4. Numerical Results 

Let us now discuss the numerical implications 
of the calculation presented in the preceding sec- 
tions. We will present the dependence of 773 on 
its various physical parameters and on the renor- 
malization scales at which particles are integrated 
out. 

773 depends on the scales \Xtw-, A<& and /t c - Fur- 
ther it is a function of the masses m tl m c and 
of the QCD scale parameter Aqcd- To establish 
a starting point let us pick a basic set of input 
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parameters 



= 0.31 GeV, 



m c (m c ) = n c 

1.3 GeV, Ah','.,, = 0.15 GeV, 

fi b = 4.8 GeV, M w = 80 GeV, 

Htw = 130 GeV, m t (m t ) = 167GeV. 



aLO 

QCD 



(29) 



In the following Aqcd is always understood to be 
defined with respect to four active flavours, the 
corresponding quantities in effective three- and 
five flavour-theories are obtained by imposing con- 
tinuity on the coupling a s at fi c and fib^\ 



The value for 773 corresponding to the set (29) 
reads: 



nf° = 0.365, 



% * NLO = 0.467. 



(30) 



Hence the NLO calculation has enhanced 773 by 
27%. From the difference of 0. 102 between the two 
values in ( |30| ) 0.022 originates from the change 
from the LO to the NLO running a s . The small- 
ness of this contribution is caused by the adjust- 
ment of Aq^, d to fit the NLO running coupling. 
The explicit 0{a s ) corrections from the NLO mix- 
ing and matching contribute 0.080. 

Let us further quantify the influence of the pen- 
guin operators Q3,...,6 : If one neglects them com- 
pletely, one obtains 7^ NLO ' np = 0.472 with the 
set in (p9|), i.e. their contribution is of the order 
of 1%. 

4.1. Scale Dependence of 773 

Ideally 773 should not depend on the factoriza- 
tion scales Htw-i ^b, fJ-tw- Yet due to the trunca- 
tion of the perturbation series such a dependence 
shows up. It may serve as an estimate of the the- 
oretical error of the calculation. 

It turns out that the dependence of 773 on fib is 
extremely mild. This is due to the fact that no 
diagrams containing internal bottom quarks con- 
tribute to the |AS|=2 process in order a s . The 
only places where fib enters are a) the running of 
a s , b) the NLO matching matrices and anomalous 
dimensions of the | AS| = 1 penguin operators. Nu- 
merically one finds that 7/3 is shifted by 0.1-0.2% if 
one choses the extreme values of fib = fi c — 0(m c ) 
or ii b = fitw = 0(m w ,m t ). 



7 ) Threshold corrections appearing for {i q ^ m q 
merically negligible. 



Now let us turn to the more important cases, 
fitw and fi c - First consider the variation of 77? 
with respect to fitw, which is displayed in Fig. |2| 

Since at fitw the top quark and the W-boson 



0.46 
0.44 
0.42 
^3 0.40 
0.38 
0.36 
0.34 
0.32 
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LO 



90 



110 130 150 170 190 
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Figure 2. The variation of 773 in LO and NLO 
with respect to the scale fitw, at which the initial 
condition is defined. The other input parameters 
are given in 

are integrated out simultaneously, it is natural 
to choose the interval Mw < fitw < "*( for the 
analysis. In the LO result for 773 we find a sizeable 
scale dependence of 12%. It is almost totally re- 
moved in the NLO, where we obtain a variation of 
less than 3% in this interval. This shows that it is 
very accurate to integrate out the two heavy par- 
ticles simultaneously. The strong improvement in 
the NLO is due to the smallness of In xt ■ 

The situation is not so nice in the case of the 
variation of fi c , which is displayed in Fig. ||. We 
have intentionally extended the range for fi c to 
the unphysical low value of 0.7 GeV to visualize 
the breakdown of perturbation theory. Varying 
fi c within the interval 1.1 GeV < fi c < 1.6 GeV 
yields 

0.33 < ?73 LO < 0.40, 0.43 < 7/* NLO < 0.50. (31) 
This corresponds to a reduction of the scale de- 
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Figure 3. The variation of 773 in LO and NLO 
with respect to the scale [i c . The range for the 
latter is taken unphysically large to visualize the 
breakdown of perturbation theory. The other in- 
put parameters are given in (f2f|). 

pendence from 20% to 14%. One reason for the 
poor improvement is the fact that the NLO run- 
ning of the mass is stronger than the LO one. 

4.2. Dependence of 7/3 on Physical Quanti- 
ties 

Let us now investigate the dependence of 773 on 
the physical parameters. From the smallness of 
the coefficient C7 at the initial scale one expects 773 
to be almost independent of m\ = mt(m t ). This 
statement is confirmed numerically, allowing to 
treat 773 as TOt-independent in phenomenological 
analyses. 

The LO result for 773 depends on to* = m c (m c ) 
sizeably. Yet this dependence is washed out nearly 
completely if one looks at the NLO 77*, see Fig. [I| 

We close this section by a look at the depen- 
dence of 773 on Aqcd j which is plotted in Fig. ||. 
It also turns out to be very moderate. 



Figure 4. The dependence of 77J on m c (m c ) in LO 
and NLO. The other input parameters are given 



5. Conclusions 

We have calculated the QCD short distance co- 
efficient 773 of the low energy | AS| = 2 hamiltonian 
in the next-to-leading order (NLO) of renormal- 
ization group improved perturbation theory. It 
reads 



n 47 +003 

u - 4 ' -0.04 1 



0.37. 



(32) 



The coefficient is scheme independent except that 
it depends on the definition of the quark masses 
in H |AS K The result in (|§) corresponds to MS- 
masses m c (m c ) and rnt(mt) as indicated by the 
superscript 

The result has passed several checks: 

i) The NLO anomalous dimension tensor j±j.7 
(21b) has been found independent of the in- 
frared structure of the two-loop diagrams. 

ii) We have kept the gluon gauge parameter 
£ arbitrary. It has vanished from 7±j,7 af- 
ter adding the contributions of the diagrams 
with their correct combinatorial weight. 
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hi) ln(m c //x)/e-terms have disappeared from 
the sum of two-loop diagrams and countert- 
erm diagrams. 



iv) The dependences of the final result for 773 on 
the matching scales HtWi Mb an d /J c cancel 
to order ot s . Numerically the dependence 
has decreased. 



v) If one expands the final result in powers 
of ot s , one recovers the terms proportional 
to a° In 1 x c , a® In x c , a\ In 2 x c , a\ In 1 x c of 
the result without RG improvement. 



The initial condition for C7 in (|l^) as well 
as the anomalous dimension tensor 7±j,7 
in ( 21b ) depend on the definition of the 
evanescent operators (0). We have checked 
that this dependence is in accordance with 
the theorems of || , so that the final result is 
independent of the choice of the evanescent 
operators. 



